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We present a study of fusion cross sections using a new generation Time-Dependent Hartree- 
Fock (TDHF) code which contains no approximations regarding collision geometry and uses the full 
Skyrme interaction, including all of the time-odd terms. In addition, the code uses the Basis-Spline 
collocation method for improved numerical accuracy. A comparative study of fusion cross sections 
for ^^O -i-i'''^* O is made with the older TDHF results and experiments. We present results using 
the modern Skyrme forces and discuss the influence of the new terms present in the interaction. 



PACS numbers: 21.60.-n,21.60. Jz 



O 
O 

(N 



> 
oo 
cn 

o 
m 
o 
\o 
o 



X 



I. INTRODUCTION 

With the increasing availability of radioactive ion- 
beams Q the study of structure and reactions of exotic 
nuclei are now possible The microscopic descrip- 

tion of such nuclei will lead to a better understanding of 
the interplay among the strong, Coulomb, and the weak 
interactions as well as the enhanced correlations present 
in these many-body systems. This has lead to a con- 
siderable theoretical effort to perform nuclear structure 
calculations with ever increasing accuracy and extensive 
investigations of the nuclear effective interaction Q . 

From a theoretical point of view, these highly complex 
many-body systems are often described in macroscopic 
terms. This has been particularly true in the case of 
non-relativistic heavy-ion collisions For example, the 
time evolution of the nuclear surface and the correspond- 
ing geometrical shape provides a very useful parameter to 
help organize experimental data. Using this approach nu- 
merous evolutionary models have been developed to ex- 
plain particular aspects of the experimental data 0,IE0|- 
These methods provide a useful and productive means 
for quantifying multitudinous reaction data. In practice, 
they require a quantitative understanding of the data as 
well as a clear physical picture of the important aspects 
of the reaction dynamics. The depiction of the collision 
must be given at the onset, including the choice of coordi- 
nates which govern the evolution of the reaction. Guess- 
ing the correct degrees of freedom is extremely hard, 
without a full understanding of the dynamics, and can 
easily lead to misbegotten results. More importantly, it 
is often not possible to connect these macroscopic classi- 
cal parameters, describing nuclear matter under extreme 
excitation and rearrangement, with the more fundamen- 
tal properties of the nuclear force. Ultimately, these dif- 
ficulties can only be overcome with a fully microscopic 
theory of the collision dynamics. 

In this paper, we utilize the time-dependent Hartree- 
Fock (TDHF) method. It is generally acknowledged that 
the TDHF method provides a useful foundation for a fully 
microscopic many-body theory of low-energy heavy-ion 
reactions 10, 11,, J^]. The TDHF method is most widely 
known in nuclear physics in the small amplitude domain, 
where it provides a useful description of collective states 



[TsL 13 [l5| . and is based on the mean-field formalism 
which has been a relatively successful approximation to 
the nuclear many-body problem for reproducing the prin- 
cipal properties of stable nuclei throughout the periodic 
table. During 1970's and 1980's the TDHF theory has 
been widely used in the study of fusion excitation func- 
tions, fission, deep-inelastic scattering of he avy mass sys- 
tems, and nuclear molecular resonances [lOl Hi lig, while 
providing a natural foundation for many other studies. 
An account of some of the previous TDHF applications 
can be found in Refs. [lOlll]]|. 

In the next section we will summarize some theoretical 
aspects of TDHF theory and give an account of earlier 
calculations as it is relevant to this work. In Section UTTI 
we present new TDHF fusion calculations and compare 
them to older results and, when available, experiments. 



II. THEORETICAL DETAILS 

Despite its wide usage, it has been difhcult to assess the 
reliability of the TDHF calculations due to an occasional 
imperfect or even incorrect reproduction of experimental 
behavior. This has naturally lead to consider various ex- 
tensions to the theory, particularly the inclusion of the 
two-body collisions However, there are im- 

portant components of the basic theory which have not 
yet been fully implemented, and the viability of the anal- 
ysis depends on the overall accuracy of the TDHF calcu- 
lations. The assumptions and approximations that may 
impact the results of the TDHF calculations can be cate- 
gorized as: (a) Symmetry assumptions about the collision 
dynamics, (b) symmetry assumptions used for the nu- 
clear force, (c) accuracy of the numerical implementation. 
Approximations of any type limit the number of degrees 
of freedom accessible during a collision, and hence the 
nature and degree of dissipation [23, |2l|, |22, UM ■ The un- 
derstanding of the dissipative mechanisms in the TDHF 
theory is vital for establishing the region of validity of the 
mean-field approximation and providing estimates for the 
importance of the mean-field effects at higher energies. In 
TDHF, the dissipation of the translational kinetic energy 
of the two ions is due to the collisions of single particle 
states with the walls of the time-dependent Hartree-Fock 
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potential. This leads to the randomization of the mo- 
tion characterized by the distribution of energy among all 
possible degrees of freedom of the system. The complete 
equilibration of the translational kinetic energy among all 
possible degrees of freedom is commonly accepted as be- 
ing the definition of fusion whereas the incomplete equi- 
libration results in inelastic collisions. 



A. TDHF Collision 

In TDHF, the initial nuclei are calculated using the 
static Hartree-Fock (HF) theory. The resulting Slater 
determinants for each nucleus comprise the larger Slater 
determinant describing the colliding system during the 
TDHF evolution, as depicted in Fig. Nuclei are as- 




FIG. 1: Schematic illustration of the initial and final many- 
body states. The initial state is block diagonal whereas the 
final state is a full Slater determinant. 

sumed to move on a pure Coulomb trajectory until the 
initial separation between the nuclear centers used in 
TDHF evolution. Using the Coulomb trajectory we com- 
pute the relative kinetic energy at this separation and the 
associated translational momenta for each nucleus. The 
nuclei are than boosted by multiplying the HF states with 

exp(ikj • R)$j , (1) 

where $j is the HF state for nucleus j and R is the 
corresponding center of mass coordinate 

R^i-Er.. (2) 

J i=l 

The Galilean invariance of the TDHF equations (dis- 
cussed below) assures the evolution of the system with- 
out spreading and the conservation of the total energy 
for the system. In TDHF, the many-body state remains 
a Slater determinant at all times. The final state is a 
filled determinant, even in the case of two well separated 
fragments. This phenomenon is commonly known as the 
"cross-channel coupling" and indicates that it is not pos- 
sible to identify the well separated fragments as distinct 
nuclei since each single particle state will have compo- 
nents distributed everywhere in the numerical box. In 
this sense it is only possible to extract inclusive (aver- 
aged over all states) information from these calculations. 

Approximations used in collision geometry include the 
assumption of an axially symmetric geometry used in 



earlier TDHF calculations ^J. In addition, reflection 
symmetry with respect to a fixed reaction plane and 
z-parity symmetry for identical systems have also been 
used. For axially symmetric calculations, non-central col- 
lisions were studied using the so called "rotating frame 
approximation" 0|. During the past decade some of 
these assumptions, specially the axial symmetry con- 
straint have been relaxed |2^ . A limited number 
of comparisons of axially symmetric TDHF calculations 
with the corresponding three dimensional calculations 
are available [25, 27, 28, 29] . The three dimensional cal- 
culations show more dissipation as anticipated. 



B. Effective interaction 

Almost all TDHF calculations have been done using 
the Skyrme interaction. A variety of calculations have 
shown that the TDHF results are very sensitive to the 
different parametrization of the Skyrme force '2Q| f2ll |2^ 
|2^ . Fusion behavior is especially sensitive to the effective 
interaction l20j . Some of the assumptions made in ear- 
lier calculations included neglecting the spin-orbit force 
and assuming spin saturation, neglect of pairing and the 
use of the "filling approximation" for the occupancy of 
the last partially filled shell, and the time-reversal in- 
variance of the single particle Hamiltonian. Most of the 
earlier TDHF calculations also replaced some of the nu- 
merically difficult terms in the Skyrme interaction with 
a finite-range Yukawa form |3Clj| , without a new fit to the 
nuclear properties. Previously, we have shown that the 
inclusion of the spin-orbit interaction lead to enough ad- 
ditional dissipation to resolve the well known "fusion win- 
dow anomaly" (a non-zero lower orbital angular momen- 
tum limit for fusion) [Hillll. Most of the new gener- 
ation TDHF programs do include at least the traditional 
spin-orbit interaction. However, it is well known '31] that 
the Skyrme energy density functional also contains terms 
which depend on the spin density, s, spin kinetic energy 
density, T, and the full spin-current pseudotensor, J, as 

E = J dVH(p,r,j,s,T, J;r) . (3) 

The time-odd terms (j, s, T) vanish for static calculations 
of even-even nuclei, while they are present for odd mass 
nuclei, in cranking calculations, as well as in TDHF. The 
spin-current pseudotensor, J, is time-even and does not 
vanish for static calculations of even-even nuclei. How- 
ever, this terms has not been commonly included in its 
full extent in the fitting of the Skyrme parameters due to 
its numerical complexity (the spin-orbit density J is the 
antisymmetric part of this pseudotensor, and has been 
included). The inclusions of these terms modifies the 
Skyrme energy density functional as. 
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where T-Lq is the Skyrme energy density functional used 
in earlier calculations, with the exception of the spin- 
orbit term containing the density J. The Skyrme en- 
ergy density functional does remain time-reversal invari- 
ant as all the time-odd terms enter in quadratic form or 
as linear bi-products. These terms, while required for 
TDHF to maintain the Galilean invariance of the colli- 
sion process |22jj have not been included in TDHF cal- 
culations because of numerical difficulty. Recently, due 
to renewed efforts towards an improved Skyrme interac- 
tion for static nuclear properties, a number of investi- 
gations have focused on identifying the importance and 
impact of these time-odd terms |^ [s^, H^- It is clear 
that they can no longer be neglected in TDHF calcula- 
tions, at least for preserving Galilean invariance and as- 
suring that TDHF calculations are truly based on the 
same static effective interaction, since the most mod- 
ern parametrization of the Skyrme force include such 
terms Finally, the pairing force has sometimes been 
included in TDHF calculations as approximated by BCS 
type pairing, where BCS equations are solved for the 
calculation of the initial static nuclei and the occupa- 
tion numbers are kept frozen during the time-evolution. 
It has previously been argued that, due to the exten- 
sive continuum coupling and internal excitations dur- 
ing the time-evolution, the effects of pairing will quickly 
wash away p^ . whereas other calculations have shown 
stronger pairing correlations • There is also the ques- 
tion of handling pairing, which is inherently related to 
time-reversal invariance, and the time-reversal breaking 
terms at the single-particle level for TDHF calculations. 
The study of the importance of pairing interactions dur- 
ing the collision process is still an open question and can 
only be properly answered by performing time-dependent 
Hartree-Fock Bogoliubov (TDHFB) calculations Fi- 
nally, most Skyrme parametrizations include a one-body 
center-of-mass correction term, which is not included in 
generating the initial static solutions for the TDHF evo- 
lution. 



C. Numerical approximations 

From the numerical standpoint, new techniques have 
been developed to handle the solution of the Hartree- 
Fock equations on a space-time lattice. Equations of mo- 
tion are obtained via the variation of the lattice repre- 
sentations of the constants of motion, such as the total 
energy 0jll3- ^^i^ approach, finite lattice equations 
which exactly preserve the constants of motion emerge 
from the theory in a systematic way. Most of the earlier 
numerical calculations have employed low order finite- 
difference discretization techniques where the resulting 
numerical accuracies limited the studies to the gross fea- 
tures of the reaction. With modern supercomputers it 
has become feasible to carry out more extensive nuclear 
structure and reaction studies employing higher-order 
discretization techniques, such as fifth and seventh order 
finite-difference. Over the last decade we have developed 
a more modern and advanced technique by discretization 
of the energy density functional on a basis-spline colloca- 
tion lattice, which provides a highly accurate alternative 
to the finite-difference method |37ll38| . 



III. FUSION CROSS-SECTIONS 

Heavy-ion fusion reactions are a sensitive probe of the 
size, shape, and structure of atomic nuclei as well as 
the collision dynamics. Fusion studies using neutron- 
rich nuclei are becoming increasingly available. In re- 
cent experiments with heavy neutron-rich ^'^'^Sn beams 
on ^'^Ni 0, enhanced fusion-evaporation cross sections 
have been observed. Another experimental frontier is the 
synthesis of superheavy nuclei in cold and hot fusion re- 
actions n41, 4^ ,;43. ^] . Some phenomenological mod- 
els predict that the fusion cross sections depend on the 
heavy-ion interaction potential and on the nuclear form 
factors in the vicinity of the Coulomb barrier 39] . The 
more recent theoretical approaches for calculating heavy- 
ion fusion cross sections may be grouped into three major 
categories: a) barrier penetration models [391, [isl [io. l47j . 

b) coupled-channels calculations [isl l49l ISOllsiLIS^ . and 

c) microscopic many-body approaches such as the TDHF 
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method ES m 111 m . 

In fusion, the relative kinetic energy in the entrance 
channel is entirely converted into internal excitations of 
a single well defined compound nucleus. In TDHF theory 
the dissipation of the relative kinetic energy into internal 
excitations is due to the collisions of the nucleons with the 
"walls" of the self-consistent mean- field potential. TDHF 
studies demonstrate that the randomization of the single- 
particle motion occurs through repeated exchange of nu- 
cleons from one nucleus into the other. Consequently, the 
equilibration of excitations is very slow and it is sensitive 
to the details of the evolution of the shape of the compos- 
ite system. This is in contrast to most classical pictures 
of nuclear fusion, which generally assume near instan- 
taneous, isotropic equilibration. Although fusion reac- 
tions occur for light, medium, and heavy-systems there 
are qualitative and quantitative differences among these 
systems. The interpretation of fusion reactions in terms 
of a semi-classical TDHF theory exhibits the best agree- 
ment with experiment for the lightest systems, since here 
fusion comprises almost the entire reaction cross section. 
Since TDHF is a semi-classical theory it is only possi- 
ble to calculate fusion cross-sections above the barrier, 
which is dynamically determined and may be different 
than the one calculated using a static two-center model. 
Historically, TDHF calculations have been shown to re- 
produce the general trends of the observed fusion data 
[Hill El. The TDHF fusion cross-section is calculated 
using the sharp-cutoff approximation |29j| 
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where ^ is the reduced mass, Ecm is the initial center of 
mass energy, and the quantity Imax denotes the maxi- 
mum orbital angular momentum for which fusion occurs. 
Previously, the above expression for fusion cross-section 
contained a non-zero lower limit for orbital angular mo- 
mentum to accommodate for central transparency ob- 
served for some systems. The so called "fusion- window 
anomaly" , which had not been experimentally seen and 
has been considered to be the breakdown of the mean- 
field approach, has been shown to disappear when the 
spin-orbit interaction was included in the TDHF calcu- 
lations [20| . 

We have carried out a number of TDHF calculations 
for the i^O-h^^O system using different parametrizations 
of the Skyrme force and compared them to earlier calcu- 
lations. The calculations were done in an unrestricted 
three-dimensional geometry using a basis-spline colloca- 
tion lattice of (—14, -t-14)^ and a lattice spacing of 1.0 fm. 
The static solutions were obtained using the gradient it- 
eration method to an energy convergence of 1 part in 10^^ 
and the time evolution used the exponential expansion of 
the infinitesimal propagator for up to 15 terms. With- 
out assuming time reversal invariance each single particle 
state is represented by a two-spinor carrying an occupa- 
tion number of n = 1. So, for a single ^^O nucleus we 
have 16 single-particle states, each having a spin-up and 



TABLE I: Calculations of the fusion threshold energy for 
^^O O using various parametrizations of the Skyrme in- 
teraction. The subscript Y indicates that the V^p terms are 
replaced by a finite-range Yukawa form for computational rea- 
sons. T = indicates no time-reversal symmetry for the in- 
teraction was assumed. 



Force 



Ethreshold (MeV) 



Comment 



Sklly 

SkMy 

SKM* 

SkM* 

Sly4 

Sly4 

Sly5 



70 
62 

56 
56 
53 
55 



Ref. [20j, 2D, only 
Ref. [20], 2D, only 
T = 0, 3D, only ^ 
T = 0, 3D, include 
T = 0, 3D 
T = 0, 3D, include 



T = 0, 3D, fitted with 



a spin-down component. The nuclei were initialized as- 
suming that they approach each other asymptotically on 
a Coulomb trajectory. The initial separation of the nuclei 
for TDHF calculations was taken to be 15 fm. Further 
numerical details and the accuracy of our calculations 
have been discussed in Ref. • 

We first examine the threshold energy for fusion, which 
is the energy above which we do not observe fusion 
but only inelastic collisions. This is done by perform- 
ing head-on (zero impact parameter) collisions for vari- 
ous parametrizations of Skyrme the interaction and com- 
pared to earlier TDHF calculations. It should be noted 
that for head-on collisions the reduction of TDHF equa- 
tions to axial-symmetry is almost exact. The results are 
tabulated in Tabled along with comments indicating the 
details of the force selection in each case. Since the in- 
clusion of the spin-orbit interaction was found to have 
a profound impact on these results ^ we do not dis- 
cuss prior results that do not include this contribution. 
The first two threshold values denote the calculations 
done using axially symmetric geometry and the Yukawa 
finite-range approximation for the Skyrme parametriza- 
tions Skll 54] and SkM* Traditionally, all TDHF 
calculations included the time-odd current j appearing in 
combination (pr — j"). The next value is the same calcu- 

TABLE II: Calculations of the fusion cross section for ^^0+^^ 
O using various parametrizations of the Skyrme interaction. 



Force 



o fusion (mb) 



Sklly [20] 
SkM;,[20j 
SkM* 
Sly5 

Experiment 



1694 
1822 
1368 
1347 
1075 



lation performed using the exact form of the SkM* inter- 
action and in three-dimension. As we see, the threshold 
energy is reduced by 8 MeV. Since, axially symmetric ge- 
ometry is almost exact for head-on collisions this differ- 
ence could be largely attributed to the incorrect surface 
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behavior generated by the Yukawa approximation and 
perhaps to substantially improved numerical accuracy. 
The next four threshold values include all of the time- 
odd terms in the Skyrme interaction. However, there 
is still an unresolved issue regarding the terms contain- 
ing the time-even pseudotensor J. This term is non-zero 
for static calculations but has not been fully included in 
most fits for the Skyrme interaction. On the other hand 
it may be necessary to maintain the Galilean invariance 
of the TDHF evolution. Repeating the same calculation 
for SkM* but including all of the terms in the Skyrme in- 
teraction results in a reduction of the threshold energy by 
another 6 MeV. Finally, we have performed calculations 
with more modern Skyrme forces, SLy4 and SLy5 |34j . 

The parametrization SLy4 does not include the con- 
tribution to the Skyrme energy density functional. The 
inclusion of this term results in a 3 MeV reduction in 
threshold energy. This is interesting because the con- 
tribution of this term to the binding energy of the ^^O 
nucleus is on the order of a few tens of keV. The last row 
of Tabled shows the result for the SLy5 parametrization, 
which includes the contribution in determining the 
force parameters. We can conclude that, despite small 
differences, most modern forces seem to yield threshold 
energies that are in agreement with each other. 
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FIG. 2: (Color online) The contribution of time-odd terms 
to the total energy is plotted as a function of time for the 
i6q ^16 Q system at fe = 6.6 fm, using the SLy5 interaction. 
On the same figure we also separately show the contribution 
arising from the time-odd term, which was also present in 
earlier TDHF calculations. 



We have also performed fusion calculations for the 
^^O O system at a center of mass energy of 34 MeV. 
The reasons for choosing this particular collision energy 
is due to the availability of older calculations and data, 
as well as increased sensitivity to the details of the nu- 
clear interaction for this relatively high energy col- 
lision. The results are tabulated in Table for various 
parametrizations of the Skyrme force. The maximum im- 
pact parameter for fusion was searched in 1.0 fm intervals 



until no fusion was observed, which was then followed by 
a more precise search in intervals of 0.05 — 0.1 fm. Max- 
imum impact parameters were found to be 6.65 fm and 
6.60 fm for SkM* and SLy5, respectively. Again, we see 
substantial improvement with the more modern Skyrme 
forces when no approximation in geometry and interac- 
tion is used. The reduction in the total fusion cross sec- 
tion of 500 mb is in the right direction but still over- 
estimates the experimental cross-section ^ by about 
25%. In order to better understand the contribution of 
the various new terms contained in the time-odd part 
of the interaction we have plotted the total energy aris- 
ing from the time-odd part of the Skyrme energy density 
functional in Fig|21 On the same figure we also sepa- 
rately show the contribution arising from the time-odd 
term, which was present in earlier TDHF calculations. 
As we see, the total contribution traces the behavior of 
the contribution from the term with a slight overall 
shift. In Fig. Owe plot the time-evolution of the contri- 
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FIG. 3: (Color online) The contribution of the time-even 
pseudotensor term to the total energy is plotted as a function 
of time for the ^^0-|-^^0 system at 6 = 6.6 fm and h = 6.7 fm, 
using the SLy5 interaction. 



bution to the total energy arising from the term. We 
have plotted two impact parameters, first for h = 6.6 fm, 
where fusion occurs, and the second at 6 = 6.7 fm, where 
there is no fusion. We observe that when the nuclei are 
far apart the contribution is small and it grows as nuclei 
approach each other. However, we see a major difference 
between the case for which there is fusion and the case 
were no fusion occurs. In the former, the contribution 
rises rapidly and reaches its maximum around the time 
of greatest overlap during the neck formation. It then re- 
mains significant during the formation of the compound 
system. In the latter case the contribution remains small 
throughout the collision process, and finally the two frag- 
ments, albeit excited, come apart and move away from 
each other. This is an interesting result since this term 
does not seem to make a major contribution to the bind- 
ing energy and does not significantly alter the parameters 
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FIG. 4: (Color online) TDHF time-evolution for the ^^O O collision at an impact parameter of b - 
fusion region, using the SLy5 interaction. The initial energy is Ecm = 43 MeV. 



7.6 fm, just above the 



of the Skyrme force when it is included in the fits |34| . 
We can conclude from this that the nuclear properties 
used in fitting the parameters of the Skyrme force is not 
triggering the physical significance of this term. 

We have repeated these calculations for the collision 
of ^^O with the neutron-rich ^^O nucleus at Ecm = 
43 MeV, which maintains the same initial velocity as 
in the ^^O O case. The maximum impact param- 
eter for fusion is found to be 7.5 fm, which results in a 
cross-section of approximately 1767 mb. This value scales 
well with the mass number and does not indicate an en- 
hanced fusion cross-section for this neutron-rich system. 
Of course, the energy is relatively high and such enhance- 
ments may be seen at lower energies or below the bar- 
rier. In Fig. ^ we show this collision for b — 7.6 fm, for 
which there is no fusion. In this deep-inelastic collision 
the final translational energy of the separating ions is 
about 20 MeV, indicating that 23 MeV was utilized for 
internal excitations. The final fragments, besides being 
highly excited, show an exchange of approximately two 
protons and a neutron to the heavy fragment. The anal- 
ysis of the various contributions arising from the terms 
in the Skyrme interaction show a similar behavior to the 
16(5 _^ 16 Q system. 



IV. SUMMARY AND OUTLOOK 

The evolution of the TDHF studies of heavy ion re- 
actions is marked by increasingly sophisticated calcula- 



tions, trying to eliminate as many of the assumed sym- 
metries as possible. This progress has closely paralleled 
the advances in computer technology. 



We have presented calculations using a new generation 
TDHF program, which makes no assumptions regarding 
the collision geometry nor the Skyrme interaction, and 
uses advanced numerical methods for improved accuracy. 
We have compared the new results with earlier TDHF 
calculations and analyzed the infiuence of the new terms 
in the effective interaction, specifically the new time-odd 
terms and the spin-current pseudotensor contribution. 
In general, unrestricted calculations and new Skyrme 
parametrizations lead to substantial improvements of fu- 
sion results. The substantially different results obtained 
by earlier parametrizations of the Skyrme force seems to 
have converged to very similar outcomes for the mod- 
ern parametrizations, a sign of major strides made in 
improving the Skyrme interaction. On the other hand, 
we find that some of the new terms make an apprecia- 
ble contribution during the dynamical evolution, while 
being absent or minimally important for the static calcu- 
lations. This suggest that improvements to the Skyrme 
parametrization are still possible by incorporating dy- 
namical features into the fitting process, along the lines 
of Refs. [s^jE^I • It seems as if the mean-field approach has 
not yet been fully exhausted, and improved TDHF calcu- 
lations may display more realistic features for heavy-ion 
collisions at low and medium energies. 
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